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FOREWORD 

Two  years  of  TIROS-N  satellite  data  are  available  for  analysis 
toward  the  development  and  testing  of  algorithms  to  classify  and 
enhance  ocean  fronts  in  the  imagery.  Unfortunately,  the  mid- 
infrared  data,  which  should  yield  good  nighttime  sea  surface 
temperatures,  are  contaminated  by  a  considerable  amount  of  multi - 
band  noise. 

The  procedure  outlined  in  this  report  will  sufficiently  reduce 
the  noise  to  accomplish  the  above  tasks. 
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C.  H.  BASSETT 
Captain,  USN 
Commanding  Officer 
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I.  INTRODUCTION 


The  TIROS-N  satellite,  put  into  service  in  1978,  was  the  initial 
operational  system  of  a  series  of  third  generation,  polar-orbiting,, 
environmental  satellites.  Of  interest  to  oceanographers  was  the 
capability  of  computing  sea  surface  temperature  from  the  mid- 
infrared  (3.55  -  3.93  pm)  and  the  far-infrared  (10.5  -  11.5  ym) 
channels  and  of  detecting  certain  surface  features  associated  with 
water  masses  and  coastal  processes  from  all  four  available  channels. 

Of  specific  interest  to  the  Navy  was  and  still  is  the  detec¬ 
tion  of  fronts  (boundaries  between  water  masses)  like  the  Gulf 
Stream.  With  two  infrared  channels  for  obtaining  surface  tempera¬ 
ture,  the  TIROS-N  data  seemed  most  appropriate  for  developing  and 
testing  an  automated  frontal  detection  system  for  fleet  applications 

One  major  shortcoming  of  this  satellite  was  the  considerable 
noise  present  in  the  mid-infrared  channel  due  to  a  malfunction  in 
the  satellite's  power  system.  The  purpose  of  this  document  is  to 
describe  a  filtering  scheme  to  eliminate  or  at  least  to  sufficiently 
reduce  the  noise  to  allow  for  any  enhancement  of  the  image  required 
to  accurately  locate  the  fronts. 

While  the  demise  of  this  satellite  in  November  1980  ends  the 
utility  of  this  scheme  for  real-time  analysis  (assuming  the  defect 
is  nonrepeatable) ,  there  are  some  practical  applications  of  archived 
data  for  the  detection  of  fronts.  For  instance,  the  testing  and 
development  of  algorithms  for  the  classification  and/or  enhancement 
of  oceanic  features  and,  specifically,  their  reliability  when 
applied  to  different  locations  and  to  images  with  varying  amounts 
and  types  of  cloud  cover  are  examples. 


II.  PROCEDURE 


The  data  processing  was  performed  on  the  Navy's  Interactive 
Digital  Satellite  Image  Processing  System  (IDSIPS)  described  by 
Pressman  and  Holyer  (1978).  Existing  software  described  by  the 
Stanford  Technology  Corp.  (1978)  was  adequate  for  the  analysis. 

The  first  set  of  data  consisted  of  a  512-by-512  pixel  (picture 
element)  image  of  the  western  North  Atlantic  containing  the  Gulf 
Stream  (fig.  1).  Applying  a  linear  intensity  enhancement  to  the 
image  provides  no  apparent  benefit  towards  interpretation  and,  in 
fact,  degrades  edge  detail  as  shown  in  figure  2.  A  pseudo-color 
enhancement  with  16  colors  (not  shown)  assigned  to  the  256  grey 
levels  of  the  image  accentuates  the  noise  to  a  degree  where  the 
edges  are  almost  undetectable. 

Figure  3  is  a  cross-section  or  profile  from  left  to  right, 
centered  in  the  image.  The  grey  level  scale  appears  on  the 
ordinate . 


The  frequencies  where  the  noise  occurs  were  located  by  examina¬ 
tion  of  figure  4  which  is  the  natural  logarithm  of  the  Fast  Fourier 

Transformed  (FFT'd)  image.  Normally,  a  finite  data  set  is  tapered 
at  the  ends  to  prevent  distortion  in  the  spectrum  due  to  the  dis¬ 
continuity  in  the  replicated  data.  However,  some  experimentation 
indicated  that  the  512-by-512  data  set  was  large  enough  to  appro¬ 
ximate  an  infinite  set  of  data  and  the  tapering  was  excluded  from 
the  scheme.  Note  that  only  half  of  the  loge  FFT  is  displayed  since 
the  missing  half  is  symmetric  and,  therefore,  redundant. 

Tht  transformed  image  exhibits  a  very  narrow  band  of  low  fre¬ 
quency  noise  located  at  .015  cycles/pixel  in  the  horizontal  and 
uniform  in  the  vertical  and  a  much  broader  band  of  high  frequency 

noise.  Figure  5  is  it  profile  of  the  loge  FFT  along  the  vertical 

axis  at  zero  cycles/pixel  on  the  horizontal  axis.  Figure  6  is  a 
profile  crossing  the  same  image  at  the  top.  In  this  figure,  the 
low  frequency  noise  is  barely  discernible  while  the  wider  band  of 
high  frequency  noise  is  clearly  evident  at  the  right  end  of  the 
prof i le . 


Cleaning  up  the  image  was  accomplished  by  applying  two  Gaussian- 
shaped  filters  to  the  two  bands  of  noise  described  above.  The 
characteristics  of  the  filters  designated  as  FILTER  1  and  FILTER  2 
are  outlined  in  table  1.  Filtering  by  the  more  sophisticated  local- 
sr.ailsti.es  method  described  by  Lee  (1980)  and  Willman  (1980)  was 
rejected  because  it  appeared  to  be  applicable  only  to  high  frequency 
no i so  . 


The  lev/  frequency  noise  was  attacked  first  with  FILTER  1  applied 
to  the  image  in  the  transform  domain.  Figure  7  is  a  top  view  of 
the  filter  and  shows  its  narrow,  elliptical  shape.  A  complementary 
view  across  the  center  of  the  filter  is  given  in  figure  8  which 
demonstrates  its  inverse  character. 

After  taking  the  Inverse  Fourier  Transform  (IFT)  of  the  filtered 
image,  the  resulting  image  was  enhanced  in  the  same  fashion  as  the 
initial  image  and  is  displayed  in  figure  9  along  with  a  profile  of 
grey  values  (fig.  10)  at  approximately  the  same  location  as  the 
prof ilo  in  figure  3.  The  lFT'd  image  without  enhancement  (not 
shown)  appeared  to  be  clear  of  all  noise,  but  figures  9  and  1 0 
indicate  the  presence  of  the  high  frequency  noise.  Note  that  the 
wave-like  pattern,  most  obvious  above  the  cloud  cover  in  the  lower 
rig, lit  corner  of  figures  1  and  2,  is  missing  in  figure  9  proving  the 
effectiveness  of  the  scheme  so  far. 

Applying;  FILTER  2,  described  in  table  1  and  shown  in  figure  11, 
reduced  the  apparent  high  frequency  noise  which  is  difficult  to 
detect  in  the  enhanced  image  in  figure  12.  For  comparison,  the 
profile  in  figure  13  crosses  the  i tnuge  at  the  same  location  as  the 
profiles  in  f’igures  3  and  LO .  The  demarcation  between  land  and 
water,  different  water  masses  (indicated  by  FRONT),  and  water  and 
clouds  (at  I  he  extreme  right  side)  is  now  obvious. 
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Table  1.  GAUSSIAN  FILTER  CHARACTERISTICS 


Parameter  Description  * 

FILTER  1 

FILTER  3 

Gain  of  filter  at  zero  frequency 

0. 

1  . 

Gain  of  filter  at  Nyquist 
frequency 

1  . 

0. 

Frequency  along  major  axis  of 

filter's  elliptical  contours  at 
which  filter  attains  Y?  its 
maximum  power 

.004  cycles, 
pixel 

.20  cycles/ 
pixel 

Frequency  along  minor  axis  of 
filter's  elliptical  contours  at 
which  filter  attains  %  its 
maximum  power 

.495  cycles/ 
pixel 

.30  cycles/ 
pixel 

Horizontal  (sample)  coordinate/ 
frequency  of  filter  center  in 
transform  domain 

10  samples/ 

.015  cycles/ 
pixel 

0  samples/ 

0  cycles/ 
pixel 

Vertical  (line)  coordinate/ 
frequency  of  filter  center  in 
transform  domain 

357  lines/ 

0  cycles/ 
pixel 

257  lines/ 

0  cycles/ 
pixel 

Angle  relative  to  horizontal 
direction  of  filter  orien¬ 
tation  in  transform  domain 

0  degrees 

0  degrees 

*  Based  on  a  frequency  scale  from  -0.5  cycles/pixel  to  +0.5 
cycles/pixel . 
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An  attempt  to  remove  the  high  frequency  noise  by  convolving 
the  image  (in  the  spatial  domain)  with  a  3-by-3  pixel  kernal  with 
all  elements  in  the  kernal  weighted  equally  was  less  successful. 
While  the  convolved  image  in  figure  14  differs  only  slightly  from 
the  twice-filtered  image,  the  profile  of  the  convolved  image 
(fig.  15)  is  less  demonstrative  in  delineating  features. 

III.  CONCLUSION 


a  scheme  for  noise  removal  in  the  3.55  -  3.93  pm  infrared 
channel  of  the  TIROS-N  satellite  is  proven  to  be  effective. 

Low  frequency  noise  is  removed  with  a  narrow,  inverse  Gaussian 
filter.  High  frequency  noise  is  removed  by  either  a  low-pass 
Gaussian  filter  or  by  convolving  the  image  in  the  spatial  domain. 
The  latter  technique  is  less  desirable  since  edges  are  smoothed 
over  somewhat  and  some  high  frequency  noise  is  still  present  to 
clutter  the  fine  detail. 

This  scheme  was  applied  to  another  noisy  image  with  similar 
results.  Therefore,  the  process  is  nonunique. 
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Profile  of  grey  levels  across  the  center  of  the  image  in  figure 
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